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Abstract. Ultra-cold alkali atoms trapped in two distinct hyperfinc states in an 
external magnetic field can mimic magnetic systems of spin 1/2 particles. We describe 
the spin-dependent effective interaction as a spin-spin interaction. As a consequence 
of the zero-range, the interaction of spin 1/2 bosons can be described as an Ising or, 
alternatively, as an XY-coupling. We calculated the spin-spin interaction parameters 
as a function of the external magnetic field in the Degenerate Internal State (DIS) 
approximation. We illustrate the advantage of the spin-spin interaction form by 
mapping the system of N spin 1/2 bosons confined by a tight trapping potential on 
that of N spin 1/2 spins coupled via an infinite range interaction. 
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1. Introduction 

The ability to mimic quantum spins and the promise for controlhng their mutual 
interactions are central ingredients of the cold atom prospects for simulating complex 
quantum systems: Quantum computation schemes rely on the entanglement of qubits 
[1], which act as spin ^ objects. Lattice spin models have provided successful 
descriptions of strongly correlated electron systems [2], usually with the parameters 
situated near a quantum phase transition so that the systems exhibits quantum critical 
behavior. A general description of quantum phase transitions has not been devised 
yet [3] and while scaling arguments have been applied [1], a general renormalization 
scheme remains an outstanding challenge. The experimental inaccessibility of the 
condensed strongly correlated electron systems certainly contributes to the challenge of 
understanding quantum phase transition physics. Cold atom experiments, which offer a 
very different accessibility and a very different set of control knobs may offer a different 
and independent laboratory to test concepts. Spin | degrees of freedom are also essential 
building blocks of quantum fluids in regimes that are novel for table-top experiments 
such as Cooper-paired fermion fluids in the BEC-BCS crossover regime [El El [7], and 
polarized fermion superfiuids (for example, see theory: [HI El [ini EH [121 ESI [HI [151 [IH] 
and experiments: [T71 [HI [191 [201 [2ll [22]-) In addition, bosons trapped in two hyperfine 
states can act as spin | particles [23]. In Bose-Einstein condensates (BEC's) of such 
particles, the formation of spin ^ magnetic domains has been observed [2lj, which is a 
phase separation transition |25j . 

The observation of spin domain formation [21], of intricate multi-component spin 
physics in BECs [26], of unstable spin dynamics in spin 1 BECs [2Z] and the engineering 
of two-component BEC vortex states [2H] and [22] highlight the cold atom ability 
to access, populate and control hyperfine spin states in cold atom traps. Other 
experiments demonstrated the cold atom ability to select specific hyperfine spin states 
in a static, external magnetic field as effective spin components: The above mentioned 
observation of BEC-phase separation, the probing of the BCS-BEC crossover regime, 
the spin polarized fermion superfiuids, and the recently observed phase transition in 
an itinerant ferromagnet-like system were all carried out with effective spin ^ particles 
that occupy a linear combination of hyperfine spin eigenstates in the presence of a static 
external magnetic field. For reasons we describe below, we refer to the effective spin as 
'pseudospin'. 

In this paper, we derive the short-range, effective pseudospin-spin interaction 
potential that describes s-wave interactions of ultra-cold atoms that occupy a 
superposition of two hyperfine states in an external magnetic field. The resulting spin- 
spin interaction has an anisotropic form, but the commutation or anticommutation 
relations that respectively characterize bosonic and fermionic systems, combined with 
the short-range (delta-function) nature of the interaction project the interactions onto 
the pseudospin triplet subspace for bosons and on the pseudospin singlet subspace for 
fermions. For bosonic atoms, we find that the inter-particle interactions can be described 
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by three terms: the first term is a short-range spin independent interaction potential; the 
second term describes the coupling of one particle's pseudospin to a short-range, spin- 
independent effective magnetic field carried by the other particle; and the third term is 
a short-range Ising spin-spin coupling. Alternatively, the spin-spin interaction can be 
cast into the form of an XF-coupling. We show how these terms can be calculated in 
the Degenerate Internal State (DIS) approximation [30]. We illustrate the advantage of 
the spin-spin interaction description by deriving the many-spin hamiltonian of N boson 
particles contained in a tightly confining trap. This system, a controllable quantum 
magnet, is a promising system to probe macroscopic quantum tunneling, realize spin 
squeezing and Heisenberg- limited interferometry [3T] . 

2. Cold atom spins 

The first question to be addressed in simulating one system by another is: what are the 
salient features of the simulated system that the simulating system needs to possess? As 
we see below, a spin operator has three components that satisfy a particular commutator 
algebra. Any operator that satisfies this algebra is a candidate to acts as a spin. We 
construct the available alkali atom spin states in an external magnetic field explicitly. By 
selecting two of these states the experimentalists create a cold atom system of effective 
spin i particles. 

2.1. General spin properties: spin and rotation 



We denote the quantum state of a spin 4 by the spinor |,^). A spin of amplitude ^, 
where ^ is an integer or half-integer number, has 2^ + 1 eigenstates of the spin 
projection operator = j with j = — ^, — ^ + 1, — 1,^. The spin operator 

is the generator of rotations so that under a rotation TZ^^a by an angle a around the 
direction of unit vector e, the spin state |^) (which is a superposition of the states 
transforms according to 



Quantum mechanically, an operator O in the Hilbert space of a single spin is 
characterized by its matrix elements {C,i\0\C,j) where is the i-th element of the 
{^^, ...}, basis. We relate the basis matrix elements in the original and in the 
rotated spinor basis IC.I). When the rotated frame is obtained by rotating around e over 
an angle —6a, 




(1) 



For an infinitesimal rotation, a = 5a ^ 1, to lowest order in 6a, 




(2) 




(3) 
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where the subscript (— ) outside the bracket indicates the usual commutator, 
consequence, the operator O transforms as 



O' = 0-i6a 



As a 



(4) 



is 



Moreover, if the O-operator is rotationally invariant, O' = O, then O, T" ■ ^ 
the necessary and sufficient condition to ensure that the (9-operator is invariant under 
a rotation around the "e^-vector. 

Classically an infinitesimal rotation transforms a general vector 17" to "^-^ {'^) = 
if + (5« T" X 17". If spin is to act as a vector under a rotation, it is to transform as 

'^^M {t)=t+6a^xt- (5) 

Choosing O = , by combining Eqs. (j4]) and (E]), and identifying the terms linear i 
6a we obtain 



m 



i~t X 



The scalar product of the above vector equality with a unit vector if gives 



X ^] 



(6) 



(7) 



and choosing the ff and e vectors as part of a Cartesian (x, y, z) reference frame, ff = x 
and e = y , we obtain 



(8) 



the usual form of the angular momentum commutator relation. From this commutator, 
it can be shown that ^x^^y are raising and lowering operators, increasing and decreasing 
the spin projection eigenvalue of by one unit. As operators evolve according to the 
Heisenberg commutator equations in the quantum evolution of the system, any operator 
satisfying the above commutator relations will evolve as a spin (or angular momentum) 
of the same Hamiltonian. 

The pseudospins we describe below are not the generators of rotation and their 
mutual interactions are not rotationally invariant. However, the pseudospins can be 
described by pseudospin operators that do satisfy an angular momentum algebra. In 
the Heisenberg picture, the quantum equations of motion shows that the dynamics is 
identical to that of real spins except that the spins and their mutual interactions are 
not rotationally invariant. 



2.2. The spin structure of alkali atoms in a magnetic field 

Alkali atoms have two spin variables: the electron spin Sg of magnitude Se = | carried by 
the s-wave valence electron and the nuclear spin i of magnitude i. The electron 'spins' in 
the magnetic field of the nucleus and the short-range part of the nuclear spin's magnetic 
field (the Fermi contact term) contributes to the energy in a first-order perturbation 
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calculation [32] • The corresponding contribution takes the form a^fi ■ Sg, which is the 
hyperfine interaction. We refer to ahf as the hyperfine energy. 

In a static, external, homogeneous magnetic field B of strength 5, B = Bz, the 
interaction of the electronic and nuclear spins with the magnetic field and with each 
other are described by the Zeeman spin Hamiltonian 

T-L = Bz ■ [fieSe - fJ^Ni] + ClhfSe ■ i, (9) 

where /ig = h/2meC, fiN = h/2mpC represent the electronic and nuclear Bohr magneton, 
respectively, with rUg for the electron mass and nip for the proton mass. As the nucleons 
are heavier and the nuclear Bohr magneton, fi^, proportionally smaller (by three orders 
of magnitude) than fie, f'N and we neglect it for now. Scaling the magnetic field strength 
in units of the hyperfine field Bhf = atf/He, b = B/Bhj, the Zeeman Hamiltonian takes 
the form 

= Sg ■ i + 6 z • Sg. (10) 

tthf 

At zero magnetic field, 6 = 0, the structure of the spin Hamiltonian eigenstates is most 
readily analyzed by introducing the total hyperfine spin operator, f = i + Sg. The 
Hamiltonian eigenstates can be chosen to be eigenstates of good (f^) and f ■ z quantum 
numbers. The eigenvalues of the operator are / (/ + 1), with / = i ± 1/2. The 
eigenvalues of the f ■ z are labeled by the eigenvalue mj. The \f,mf) states of the 
same /-quantum number are degenerate. The first term of Eq. ffTOl) can be written as 
(f2 -si- i2) /2, and when 6 = 0, E{f,mf,b) = Eo{f) and Eo{f = i + l/2) = ahfi/2 and 
Eo{f = i — 1/2) = -aiifii + l)/2, giving a hyperfine splitting a^fii + 1). This splitting 
is determined spectroscopically. For example, "^^Na has a hyperfine energy measured to 
be ahf = 4:2.5mK = 0.95GHz (which greatly exceeds cold atom trap depths ~ fiK), 
corresponding to a hyperfine magnetic field strength equal to Bhf ^ 709G. 

At finite magnetic field, / is not a good quantum number, but it is convenient 
and customary to refer to the /-value that the same state takes on in the adiabatic 
b — 7- 0-limit. On the other hand, the total spin projection, mj remains a good quantum 
number. We parametrize the corresponding state by an angle 6 as 



\f^,mf) = ±cos (^y^ 



1 1 

mi = nif - -,mg = +- 



-|- sm — 



mi = m/ + ^,mg = -iy (11) 

We restrict the range of ^-angles to the [0, vr] interval, allowing for coefficients of a 
different sign by adding ± explicitly to the cosine term. By parametrizing the state as 
in Eq. ffTTj) , the 9 is also the inclination angle of the average electron spin vector in that 
hyperfine state. The optimal value of 9 depends on the spin projection eigenvalue mj. 

Using the angular momentum algebra Eq. ([8]), the spin product expectation value 
(Sg ■ i) , of the states of Eq. (fTTl) give 



cos {9)mf ± sin (9) 
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where, ± refer to the relative sign of the electron spin up and down components, as 
introduced in f|TT]) . Taking the extreme of the expectation value of Eq. f lTOj) with respect 
to 6 in ( fTTl) gives 

""^^ ={b + mj) cos {6) ± sin {6) + ^ ~ ^/ ~ (^3) 
and by requiring dE/dO = 0, we obtain the optimal inclination angle tan (9^), 



m 



2 



tan ie^) = ± ^ \^ . . (14) 

As 6 3> 1, tan(^='') — )■ ±[i + l]/6 ±oo, corresponding to — and Q~ — )■ tt: 
the electron spin in hyperfine states align or anti-align with the external magnetic field 
when the strength of that field significantly exceeds the hyperfine field strength. For 
notational convenience we introduce the square root of the sum of the squares of the 
numerator and denominator, 



eKm^,h) = \\\? ^ 2mfb + (^^ + ^ ) • (15) 
Then the angular projections can be written as 



Substitution of Eq. flTB]) into Eq. fITU]) yields the Zeeman eigenenergies E {±,mf,b) , 

Ei±,mf,b) ^ ^ e{mf,b) _ ^^^^ 
tthf 2 

where the superscripts ± are now seen to indicate whether the hyperfine state belongs 
to the f~^ = i + 1/2 or the f^=i — 1 /2-superpositions of Eq. ( fTTl) . To see that, 
note that in the limit of vanishing magnetic field, limfe_^05 e(mj,6) — t- ? + 1/2, so 
that E{-\-,mf,b) approaches Eq {f = i -\- 1/2) = ahfi/2 while E{—,mf,b) approaches 
Eo{f = i — 1/2) = —tthf [i + 1] /2. Actually, the expression Eq. f lTTl ) describes the 'non- 
stretched' states, |m/| < i + 1/2. For the stretched states \ f = i + 1/2, mf = i -\- 1/2) 
and \ f = i + 1/2, mf = -i - 1/2), 

E(f+,mf = ±f+,b) b i 

^ ^ ^^ = ±- + -, (18) 

ahf I I 

corresponding to 6* = and 6' = vr in the linear superpositions of Eq. f lTT]) . In Fig. [H we 
show the magnetic field dependence of the spin states for the nuclear spin value i = \ 
that describes Li, "^"^Na and ^'^ Rb, three of the most used cold atom alkali atom species. 
The electron spin expectation value (Se), 

11 . V cos (6) ib + mA , . 

(/.m;|s.|/.™,) = z^ = ±zi^. (19) 

points in the direction of the external magnetic field. 
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Figure 1. Hyperfine energy levels for nuclear spin i = 3/2, the nuclear spin magnitude 
of bosonic atoms such as "^Li, "^^Na and ^"^ Rb. 

2.3. Pseudospin 

By trapping cold atoms in a specific hyperfine state in an external magnetic field and 
by coherently transferring part of the spin population into a different hyperfine state, 
either by accessing a laser-driven two-photon Raman transition or by means of another 
magnetic field that is oscillating near the energy-difference, the experimentalists prepare 
the atoms in a superposition of two hyperfine states, |/i,m/^) and \f2,mf^). If the 
hyperfine states have been selected to ensure that no other spin states couple, |/i,mj^) 
and 1/2,^/2) act as the basis of an effective spin |. We can than assign one the role 
of 'spin-up' the other that of 'spin-down' state, |/i,mj^) = |'|') and \f2,mf^) = ||). 
Interestingly, if the atoms are bosonic, i.e. if their hyperfine spin / takes on an integer 
value, the atoms become effective spin ^ bosons. 

The coupling to other hyperfine states has to be avoided because not only does the 
role of additional hyperfine states require a larger Hilbert space description, but also it 
generally causes significant particle loss. For magnetic field strengths near the hyperfine 
magnetic field, 6 ~ 1, the Zeeman levels are widely spaced with energy differences 
comparable to ahf (tens of mK), considerably larger than the typical trap depths {^K). 
Hence, spin flip collisions in which particles end up in lower energy spin states create 
particle pairs of sufficient kinetic energy to evict spin flipped particles from the trap. 
However, conservation of energy and conservation of total spin projection (mj) in binary 
atom interactions limit the spin-changing collisions that can take place in an external 
magnetic field. 

In the case of nuclear spin i = |, the experimentalist can select either the two 
lowest or the two highest energy levels of the / = 1 manifold: either |1,0) and |1, — 1) 
or |1,0) and While colhsions |1,0) + |1,0) -)■ |1,-1) + would produce 

spin states not included in the spin- up and spin-down basis, in the region 6 ~ 1, the 
energy cost of a spin flip up {mj from to outweighs the energy gain from the 
corresponding spin flip down process (m/ from to —1), as can be seen from Fig. [T]so 
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that this process is energy forbidden. For nuclear spin 1, i = 1, as in the case of ^Li, 
^^Na and ^^Rb, the experimentahst can select the two lowest energy hyperfine states, 
It) = 1), and II) = ||, — for instance. 

In the |/i,mj^) = |t) and |/2,''^/2) = \i) basis the pseudospin operator takes the 
form s = |(T, where a is the Pauli-spin vector operator, a = a^^x + cr^y + a^z, with 
ax, (Ty, and the Pauli spin matrices. Note that s is not the generator of rotations. 
Under rotation, the spin components do not transform among each other according to 
Eq. ([T]), although the s components still satisfy angular momentum commutator relations 
[sa:,Sy]_ = isz- Therefore, the time evolution of the many-spin system (governed by 
commutation relations in the Heisenberg picture) is indistinguishable from the quantum 
evolution of the magnetic many-spin system of the same parameters. The concept of 
effective spin parallels closely that of iso-spin in nuclear physics [33], except that a 
Raman transition and/or oscillating magnetic field can convert spin-down into spin- 
up particles. Likewise, as emphasized by Bloch [31], any two-state system (and its 
decoherence) can be described by an effective spin ^ system. 

To picture the role of the pseudo-spin direction, we parametrize an arbitrary, 
normalized two component spinor \C,ps) by introducing two angles 6 and 0. 

In this notation, denotes the relative phase: the difference of the complex phase of 
the 'up' amplitude and of that of the 'down' amplitude. Also, cos6' denotes the effective 
polarization: the difference between the up probability and the down probability. Then 
the expectation value of the effective spin vector (s) in terms of 6 and (p is 

(s) = - {C^psl a \^ps) = - (sin 9 cos0x + sin ^ sin 0y + cos 9z) , (21) 

which is a vector on the surface of a sphere of radius | (known as the 'Bloch sphere'). 
The inclination 6 and azimuthal </> angles of the expectation value of the pseudospin 
vector then respectively characterize the 'polarization' and the phase difference. If an 
ensemble measurement reveals N^(i) atoms in the 'up' ('down') hyperfine state, then 
{cos{e)) = {N^-N^)/{N^ + N^). 

3. Spin-spin binary atom interactions 

Cold alkali atoms are interesting building blocks for simulating magnetic systems: their 
mutual interactions preserve the overall spin projection (the sum of m/ is preserved), 
the interactions are naturally spin dependent and the strength of these interactions and 
their spin- dependence can be varied by a Feshbach resonance. We write the effective 
atom-atom interaction as a spin-spin interaction. The statistics of indistinguishability, 
combined with the short-range nature of the effective interaction gives a short-range 
Ising-like spin-spin interaction or, alternatively, a short-range XF-interaction. We 
calculate the spin-dependence in the Degenerate Internal State (DIS) approximation. 
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3.1. Effective inter-particle interaction in pseudospin language 

In cold atom experiments, the length scales relevant to the many-body physics 
description (~ fim) significantly exceed the length scale on which the atoms interact 
(which ranges from Bohr-radius to nm). As a consequence, the interactions of 
indistinguishable atoms i and j with position coordinates Xj and Xj can be described 
by an effective contact interaction potential 

V (x,, X,) = ——a6 (x, - X,) . (22) 

where a denotes the scattering length and m represents the single particle mass. In the 
low energy regime of interest in a many-body description, this interaction reproduces 
the real binary atom scattering physics (to all orders) in a first order perturbation 
calculation. The interaction potential reproduces the correct binary atom s-wave 
scattering amplitude in the Born approximation. Introducing the annihilation (creation) 
field operators tp {ip''), the particle-particle interactions are described by 



Hint = I j d^x j rfV V'nx)^nx)^(x-x)^(xO^(x) 

d^x , (23) 



1 /Anff 



2 \ m 

in the Hamiltonian operator. When the indistinguishable particles are bosons that 
occupy two possible spin states, |t) and |4-), we distinguish interactions between particles 
in like spin states, described by a scattering length (a^) if that state is the 'up' ('down') 
spin-state, and between particles in unlike spin states, described by scattering length 
a„. The Hamiltonian operator that accounts for these interactions takes the form 







H,„, - \ ( 






m 




[ d^x 


+2a„ 





j d^x ?/'|(x)?/^j(x)?/'f-(x)?/^-|-(x) 
d'x ^A|(x)V'|(x)^/;^(x)V'^(x) 

d^x ^|(x)^/;|(x)^^(x)V;t(x)') . (24) 



By virtue of Pauli exclusion principle, fermion particles occupying two spin states W) 
and IJ,) only interact via short-ranged interactions if they are in different spin states, 

H,nt = (^^^) / ^|(x)^j(x)V;^(x)^t(x) , (25) 

where ap describes the low energy fermion- fermion scattering. We obtain these 
expressions from a scattering picture and derive a pseudo spin-spin form of the particle- 
particle effective interaction potentials. 

If the effective pseudo spin projection is conserved, interacting spin \ atoms, 1 
and 2, can undergo four types of scattering events, represented in Fig. [21 |tt) — |tt), 
III), Iti) Iti) (with the same amplitude as |it) ^ lit) ) and \\,-\) |ti), 
(with the same amplitude as W\) — l^t) ) where the first arrow in the brackets indicates 
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the spin of atom 1, the second arrow shows the spin of atom 2. By construction, 
the effective interaction yields the correct transition matrix elements in the Born- 
approximation for the four types of binary atom scattering events: (ATiffai/m) for 
the mutual scattering of two 'down' particles, (iTrh'^a^/m) for the mutual scattering of 
two 'up' particles, (47r ^^OD/m) for direct scattering (unlike spin scattering without spin 
flip) and {inh'^ax/m) for exchange scattering events (unlike spin scattering with spin 
flip). 

Denoting the full coordinate of particle i, which consists of position Xj and spin Sj 
by Tj, we can replace the effective interaction potential Eq. (l22l) by 

V{v^,r2) = a5(xi-X2), (26) 



which is a spin-operator. In the basis of the up and down spins, {|tt) > \ ti) > lit) i lii)}; 
the a-operator reads 



(27) 





a-|- 








\ 







an 


ax 












an 





v 











«4 / 



in accordance with the above description. 




1 



1. 



(a) spin up + spin up 




(b) spin down + spin down 




t i 

(c) spin up + spin down 




(d) spin up + spin dow 
(spin flip occurs) 



Figure 2. Graphic illustration of the four types of binary atom scattering processes: 
(a) two spin-up particles collide, (b) two spin down particles collide, (c) one spin up 
and one spin-down particle particle scatter while the spin projection of the interacting 
particles remains the same, and (d) a spin-up and a spin-down particle exchange spin 
states in the collision. 

To connect with spin physics we write the spin dependent scattering length 
operator, a, as the sum of products of single particle spin-operators. We write Eq. 
(EZD as 

a = aoi+a^^+ (si,^ + S2,^)+a^,p [4 (si,^S2,^)]+a^ [2 (si,^S2,x + Si,yS2,y)] , (28) 
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where I denotes the unit operator. The operators in square brackets are the spin 
projection sum operator, 

/ 1 \ 





-1 



(29) 



the spin projection product operator, 

/ 1 \ 

0-100 
0-10 




Sl,zS2,z 



(30) 



and the exchange operator, 

Sl,xS2,x + Si^yS2^y 



\ 



(31) 



/ \ 

10 
10 


By comparing the off-diagonal matrix elements, it is clear that the exchange operator 
coefficient must be equal to the exchange scattering length a^- Identifying the diagonal 
matrix elements we find Oo, ^2,+; 



a-|- + aj^ — 2a 



D 



4 



(32) 



Hence, the expression 

V{ri,r2) 



m 



ctf + + 2ao 







+ 


2 



5(xi - X2 

+ [a^ + a^- 2az)] (si.zSa,^) + 2a^ {si^^S2 
translates the effective interaction potential operator Eq. fl2B|) into pseudo-spin language 



Sl,z + S2,z) 
Sl,yS2,y)} ■ (33) 



3.2. A magnetic-like interaction form of the inter-particle interactions 

In the interest of exploiting the magnetism analogy, we note that the effective spin- 
dependent interaction Eq. (l33il can be rephrased as a magnetic interaction. Specifically, 
the spin-dependent part stems from the interaction of the effective spin of one particle 
with an effective short-ranged magnetic field carried by the other particle. That 
description allows one to interpret the interaction of one particle in a many-body system 
to be caused by the interaction with an external field and with an internal magnetic 
field generated by the moments of the other particles. 
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We write the effective interaction of atoms 1 and 2, as tlie sum of a spin independent 
contact interaction, Vo and and a spin interaction, Vg, V{fi,f2) = K (xi — 3:2) + 
Vs (ri, where 

Vo (xi - X2) = (at + ai + 2aD) S (xi - X2) . (34) 

m 

The spin-dependent part we cast into the form of the energy of a spin in a magnetic 
field with unit effective Bohr-magneton (as we chose the spin to be dimensionless, the 
effective magnetic field than takes the units of energy) 

K (ri, r2) = Si ■ h2 (xi) + S2 ■ hi (X2) , (35) 

where we assume the effective h-magnetic field to be the short-range part of an effective 
moment m, 

hi (x) = mi6 (x - Xi) . (36) 
The effective m-moment is an operator, 

nij = nis o'z + msjjZ ( Sj ■ z ) + m^^j, z x ( Sj x z ) , (37) 
with moment parameters 

ms,o = (at - «4) ) 

"^^,11 = (% + H- 2aD) , 

_L = a^, (38) 

m 

that depend on the scattering lengths. 

3.3. The effect of quantum statistics on short-range spin-spin interactions 

As the pseudo spin operators are not the generators of rotations, the spin-spin interaction 
can be and generally is anisotropic. The rrio-term describes the interaction of s with 
an effective short-range magnetic moment that points in the ^-direction (the direction 
of the physical magnetic field, the actual field that splits the Zeeman-levels) and is 
independent of the pseudo-spin of the particle that carries the moment. Another source 
of anisotropy is the difference in the my and m^-moment parameters. The short-range 
nature of the interactions and the quantum statistics of the interacting particles, give 
a moment that can be chosen to align itself with the magnetic field direction or to be 
perpendicular to it. To see that, consider the interaction 

Hint = \ y d\ (^t (^) (^) fp (^) ^ (^)^ ^ (39) 

where a denotes the spin-dependent scattering length operator of Eq. (1271) . where 

(^) ^4 (^) 
^4 (^) (^) 

V (^) (^) j 



(40) 
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indicates the two-particle spinor. The short-range nature of the interaction causes 
the field operators to be evaluated at the same position, x. As pairs of annihilation 
and creation operators of the same argument, the two-particle spinor components 
obey commutator (anti-commutator) relations if the interacting particles are bosonic 
(fermionic), {x) ipj [x) = (x) ipi [x). By writing the components of the two particle 
spinor Eq. pUl) as half the sum with itself, then replacing the second term by ± its reverse 
order, we obtain, for bosons 



ip (x) ip (x) 



for bosons and for fermions 



1 ( C'^) (^) +'^-1. (^) (^) 
1 ( V'-i. (^) (^) (^) (^) 

i^i {x) (x) 



_ \ 

' V't ('g) (^) - ^-1- (^) V't (^) ' 

V J 



(41) 



(42) 



The resulting non-vanishing column matrix elements are components of the pseudo-spin 
triplet manifold in the boson-case and the pseudospin singlet state in the fermion case. 
We find that the two-particle spinor is projected onto the triplet subspace if the fields 
are bosonic and onto the singlet subspace if the fields are fermionic. 



ip (x) ip (x) 
'ip (x) 'ip (x) 



bosons 



fermions 



'ip (x) -ip (x) 
'ip (x) -ip (x) 



(43) 



where the IIt ij^s) project onto the two-particle pseudospin triplet (singlet) subspace. 
We can check that statement by direct inspection. For instance, the second component 
of the column matrices Eq. ( I^Tl) and Eq. (H2l) represents the two-particle annihilation 
field of the | t4)-state, which can be written as | ti) = ^ (^"^^^^ + ^"^^^^j ' so that 
HtI ti) = 75 ^^'''^^^ 2ind flsl tl) = -ijJii^^ii^, corresponding to the spin states of the 
second component of the right-hand sides of Eqs. (HT]) and (H2l) . 

Writing the boson matrix product of the interaction Hamiltonian out, we find that 
we can also write the effective interaction Hamiltonian as a bracket of the triplet two- 
particle spinor, 

(^)^i (^)'/'t (^) ('44') 
(x) (x) 



■p (x) -p (x) 



V 



in terms of which 

Hint 



m 



I 



d^x T {'p'^ (x) Ip'' (x) \dT\ip {x)'ip {x] 



(45) 
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where the triplet scattering length operator clt is now represented by a diagonal matrix, 
fa-. \ 



a^ 











«4 / 



(46) 



One more application of the commutator relations casts the interaction Hamiltonian 
Eq. (145|) in the form of Eq. provided we identify the unlike boson scattering length 
au with ttu = ttD + CLx- Likewise, the fermion interaction 

1 fAnh^^ 
2 



int 



m 



d^x 5 ('j/'''^ (x) 'j/^''^ (x) lasl"?/^ (x) -j/^ (x) 



(47) 



with two-particle singlet spin component 

_ V't (^) (^) - V'i (^) V't (^) 



i/j (x) ip (x) 
and singlet scattering length 

ds = O'D — cix , 
reduces to the form of Eq 



V2 



(48) 



(49) 



if we identify ap with ap — ap) — ttx ■ 
As an interesting consequence of the triplet projection caused by boson statistics, 
the spin-spin interaction can be written in two equivalent forms. As si ■ S2 = 1/4 in a 
triplet state, we can either replace Si^zS2,z = 1/4 — {si^xS2,x + Si^yS2,y) or, alternatively, 
Si,xS2,x + ^i,y^2,y by Si^xS2,x + ^i,yS2,y = 1/4 — Si^zS2,z- As a cousequeucc of the second 
replacement, the effective spin-spin interaction takes the form of a short-range Ising-like 
interaction, the second replacement gives a short-range XY spin-spin interaction. The 
resulting boson spin-spin interactions 

A7rh^\ ^ ,^ {\a^ + a<^ + 2au 



Vi (ri,f2 



+ 



m 
a^ 



6 (xi - X2) 



(si,2 + S2,z) + K + 0-1- 2an] (si^^sa,. 



(50) 



and 



VxY (n,^2) 



S (xi - X2) 



+ 



[Sl,z + S2,z 



m 

- [at + a; - 2aJ (si,:i.S2,x + si^yS2,y)} (51) 

are equivalent to each other and to the more conventional expression of Eq. (IMl) . One 
advantage of the XF-form, Eq. (15T1) . Vxy (ri,r2), is that the unlike boson scattering 
length ttu only occurs in the spin-spin term. As a consequence a mixed spin channel 
Feshbach resonance will vary only the Xyspin-spin coupling. In general, a Feshbach 
resonance occurs when the incident particle channel becomes degenerate with the quasi- 
bound state of another collision channel. Hence a particular resonance will either vary 
a^, or ai or a^. A mixed spin channel resonance of magnetic field width AB around 
magnetic field strength Bres varies the unlike scattering length a„ as 



AB 



Pseudospin and spin-spin interactions in ultra-cold alkali atoms 



15 



where a„ is the background scattering length that varies slowly with magnetic field 
(on the magnetic field scale of B^f. The Feshbach variation thereby adds an effective 
interaction potential 

AVxY {ri,r2) = \— — — ^ I I — — — j o {xi - X2) {si^^S2,x + Si,yS2,y) {53} 

to the above XY form of the spin-dependent particle-particle interaction Eq. (IST]) . 
In terms of the magnetic form of the inter-particle interactions, 

Vi(x,Y) = ^o,/(x,y) (^1 - X2) + Vsj(x,Y) (n, ^2) , (54) 
where VoJ(xy) = [47r/i^/m]5 (xi — X2) a'o,i{x,Y) in which 

+ ai + 2au 



4 



ao,XY = — ^ — , (55) 

denote the Ising and XY expressions of the spin independent scattering length. The spin 
interactions take the usual form K,/(jfy) (^^1,^2) = ^1 ■ ^2,i{xy) i^i) + S2 ■ hi /(xY) (3:2) 
with hij(xY) (x) = 6{x - Xi) mij(xy) and 

mi,7 = m^z + (m|j - m^) z (Sj ■ z) , 

™-i,XY = TTlo'L + — my) Z X (s^ X z) , (56) 

with moment parameters (mo, mj,) defined in Eq. ( 138|) . 

3.4- Spin dependence of alkali atom interactions and the degenerate internal state 
approximation 

The spin-dependence of the alkali-atom interactions stems from the exchange of valence 
electrons. As two alkali-nuclei approach each other to nanometer and sub-nanometer 
distance, their valence electrons, now encircling the two closely-spaced nuclei, become 
strongly correlated. If these electrons are indistinguishable, i.e., if their spins align in 
a triplet sate, the likelihood of finding them simultaneously in each others vicinity is 
reduced by virtue of the Pauli principle. Pauli exclusion then reduces the Coulomb 
energy shift in the inter-atomic potential. In contrast, spin singlet electrons can 
approach each other more closely, shifting the inter-atomic potential upward. Hence, 
the triplet potential Vt is generally deeper than the singlet potential Vs- The overall 
interaction of atoms can be expressed by an inter-atomic potential operator 

V (n, V2) = (|X1 - X2I) fle.T + (|X1 - X2I) fle,5, (57) 

where the r represent both the spatial coordinates x, and spin, and IIcT , ^e,s denote the 
projection operators for the electron triplet S = 1 and singlet 5 = states. When acting 
upon a triplet state, the square of the total electron spin operator, Sg = | + 2Se,i ■ Se^2 
yields an eigenvalue of 2 (i.e. S'(S' + 1) with S = 1). Acting upon a singlet state, the 
same operator gives zero, so that the triplet projection operator takes the form 

He.T = Y = 4 + ^e.i ■ ^'=.2. (58) 
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Since He t + ^e,s = the singlet projection operator is equal to 

fle,S = 1 - fie,T = ^ - Se,l ■ 8^,2 • (59) 

With Eqs. ( 158|) and (!59|) . the interatomic interaction potential operator Eq. (1571) reads 
V (ri, r2) = ^ [3Ft (|xi - X2I) + Vs (|xi - X2I)] 

+ [Vt (|Xl - X2I) - (|xi - X2I)] Se,i ■ Se,2 . (60) 

The depth of the Vs and Vr potentials (~ electron Volt - lO'^K) greatly exceeds the 
energy of the Zeeman spin interactions ^ ahf, tens of mK). Inside the potential well, 
r < ri, the hyperfine interaction can be neglected whereas in the outer region, r > ri, the 
hyperfine interaction determines the spin state of the collision channel. The exchange 
interaction, {Vt — Vs) s^^i ■ Se^2 also falls off rapidly (exponentially) in the outer region. 
One can then calculate the scattering wavefunction while omitting the Zeeman terms in 
r < Ti and treating the exchange interaction in the region r > ri as a perturbation term. 
In lowest order perturbation the T-matrix should then have a contribution proportional 
to s*e,i-s'e,2 evaluated for the initial and outgoing channels. The spin- independent and the 
Se,i ■ Se,2-parts of the T-matrix should reproduce the correct triplet and singlet scattering 
lengths as and ar in the limit of vanishing magnetic field and hyperfine energy. In this 
approximation, the low energy i, j k, I transition matrix element then takes the form, 

Ti,j-k,i ~ {i,j\a + a-Se,i ■ Se,2\k,s} , (61) 



(62) 



m 

where a denotes the scattering length averaged over the singlet and triplet states, 
'Sax + as^ 

T 

and where a_ represents the difference scattering length, 

a^ = ax — as ■ (63) 

This approximation is called the Degenerate Internal State (DIS) approximation. While 
it was primarily designed for calculating the two-body loss rate of atoms occupying 
specific hyperfine states in an external magnetic field. In that case, if the magnetic field 
is comparable to the hyperfine field, the above formula does not work very well as the 
wavefunction of the outgoing channel is not very well approximated by the zero energy 
wave function. We are not considering lossy channels, which will not be there if the 
hyperfine states are chosen carefully as described in the previous section. We are only 
considering the cases i = k,j = 1 or i = l,j = k with i,j =t,i, for which the incident 
and final channel wavefunctions have the same low energy value. Even in that case, the 
approximation is not always satisfied, particularly for atoms that have a naturally large 
scattering length at zero magnetic field, such as "^Li |35] and ^^^Cs [36] and for higher 
magnetic field values. We expect that the treatment, may, however, yield a reasonable 
approximation for magnetic field values near such low magnetic field resonances as the 
9G resonance observed in ^"^Rb [37] . 
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In the DIS approximation, we can calculate the above defined scattering lengths 
explicitly 

=a + a_(tt |se,i-Se,2| tt), 
=a + a- (i| |se,i ■ Se,2| ii) , 

Ox =a-(ti|Se,l-Se,2|it)- (64) 

in terms of the single electron spin matrix elements Se,i- = (t |Se| t); ^e,i = (i |Se| i), 

(tt |Se,l ■Se,2| tt) = Se,f Se,t, (65) 
(11 |Se,l ■ Se,2| ii) = Sg,; ■ Sg,;, (66) 
(ti ■ Se,2| ti) = Se,t " Se,i- (67) 

The exchange spin matrix element involves spin-fiip matrix elements of the type 

(t |Se| i) = Se,t4, SO that 

(ti |Se,l ■ Se,2| it) = ^e,U " Se,4t • (68) 

The spin-spin interaction parameters, 

















= a_ 




m 




27r^2 




= a_ 




m 






ms,± 


= a_ 




m 



(Se,t + ^e,i) 



e,i) 



(69) 



then depend on the single electron spin-fiip and spin matrix elements. 

As the electron spin expectation values depend on the external magnetic field, the 
interaction parameters do as well. To express the dependence explicitly, we cast the 
expressions in parametric form, choosing the electron spin inclination angles 6^, 9^ of 
the 'up' and 'down' hyperfine states of Eq. (fT9l) as variable. For notational convenience 
we introduce external magnetic field-dependent spin factors Sm,n{b) with m and n equal 
to +1 or —1, 



Sm,n {h) = (Se,t + mSe^l) " (^e.t + nSe,l] 



COS 



(^t) + 



mcos 



(%) 



COS (9^) + n COS {9^) 



(70) 



■ s 



In addition, we introduce the exchange spin factor, (b), with (6) = ^s^^^i 
To determine its value we make assumptions about how the |t) and |i) are chosen: We 
assume that their respective m/-values differ by one unit (if not, Se^-^ = s^^^^ = 0) and 
we choose the |t) to have the highest nif. We also assume that both states are chosen 
among the states with Zeeman energy-curves that slope down at high magnetic fields 
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Figure 3. The graph shows the magnetic field dependence {b = B/Bhf) of tlie moment 
parameters as calculated in the DIS approximation for i ~ 3/2, |t) = |/ = 1,™/ = 0) 
and 1/ = l,m/ = —1). The full line shows mo/ (Trfi-^/m) a_, whereas the dashed 
line shows 2 (m|| — m±) / (Trfi^/rnj a_, where a_ represents the difference of the zero 
magnetic field triplet (ot) and singlet (as) scattering lengths, a_ = ay — as- The 
dashed lines show the cosine of the electron spin projection angle of the 'up' and 
'down' hyperfine states. The higher lying dashed line plots cos(6'4.), the lower dashed 
line plots cos (6^). It is near the 6 = 1 at which the down state electron spin changes its 
direction from parallel to antiparellel to the external magnetic field that TOq changes 
sign and |m|| — mj_\ is maximized. 



(either the / -states or the stretched electron spin-down state, \f = f~^,mf = —f'^)). 
In that case the |t)-state has the lower Zeeman energy. With this convention, we find 



(b) 



1 + cos (6't) 



cos (%) 



(71) 



where the projection of the inclination angle varies with the external magnetic field as in 
Eq. ( IT^ . In terms of the iS-spin factors, the effective pseudospin interaction parameters 
take on simple forms 



a 



— S 

m 



m 

( 

m 



-ib) 



SAb) 



(72) 



Note that all three moment parameters are proportional to the difference scattering 
length a_. 

Figure [3] plots the magnetic field dependence {b = B/Bhf as defined above Eq. 
( ITOj) ) of the relevant moment-parameters as calculated in the DIS approximation for the 
special case, i = 3/2, |t) = \ f = ^,fnff = 0) and \ f = l,mf = —1) (i.e., the second 
and third lowest Zeeman energy levels of Fig. [1]). The full line plots m^, whereas the 
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dash-dotted line plots 2(m|| — both in units of {-Kfi^a^/m). For reference, note that 
the Ising form of the spin ^ boson particle-particle interaction takes the form 

V (fi, f2) = 6{xi- X2) [Ao + moSi_2S2,2 + 2 {m\\ - m^) 51,252,^] , (73) 

where Aq denotes the spin-independent interaction strength in the Ising form, Aq = 
(irh'^/m) [a^ + ai-\- 2a„]. Note that nio changes sign near b = 0.535, so that there exists 
an external magnetic field strength at which the spin-independent effective short-range 
magnetic field carried by the interacting particles can be made to vanish (without having 
to take recourse to a Feshbach resonance). The Ising spin-spin interaction coefficient 
does not change sign but it's magnitude is maximized at a magnetic field b ~ 1.75. The 
precise values of the magnetic fields at which the interaction parameters exhibit this 
behavior may be different, but we expect the DIS-approximation to give the correct 
qualitative behavior even if the DIS-approximation is not expected to be accurate at 
higher magnetic field values. 

4. A controllable N boson quantum magnet 

The spin-spin forms of the effective inter-particle interactions reveal the analogy 
with magnetic systems. As an illustration we consider a specific system that 
promises a particularly powerful and interesting simulation of a quantum magnet: N 
indistinguishable bosons occupying two hyperfine spin states confined by a tight 

spatial potential (which could be a single well of an optical lattice). We assume that 
both I 'I') and ||) experience the same trapping potential (x) = (x) = Vr (x) of 
single-particle ground state xt (x) and single particle ground state energy ct- When 
ct exceeds all other energy-per-particle values and the system relaxed to its motional 
ground state, all of the N bosons occupy the Xr-orbital and the spatial degrees of 
freedom are 'frozen', allowing only spin dynamics. We also assume that the |t) and \\.) 
experience a coherent two-photon Raman coupling which can be effected by pulses of 
near- resonant lasers (or by means of an oscillating magnetic field). The resonant Raman 
coupling also introduces a detuning e which acts as an effective energy difference between 
It) and \l) and the Raman coupling is described by a term 

N 

E [-^^ w (1^). ^ (t i+i t), . (ii) + ^ (it), . (t i-i i), , (ii) 

TV TV 

= -2Er (t) ^ s,, + e J2 Hz, (74) 
i=i i=i 

in the Hamiltonian. In the above expression, denotes the Rabi-coupling energy 
which varies in time if the coupling is caused by a pulse. 

The premise of a trapping potential sufficiently tight to freeze out the spatial degrees 
of freedom translates into an A^-particle wavefunction of the type 

^' (fi,f2, ....,fAr) = XT (5'l) XT (^2) •••XT (^TV) | 5^ (Si, S2, ...Sat)) , (75) 
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where jS*^) denotes the spin state of the A^-boson system. By virtue of permutation 
symmetry - the full wavefunction has to be even under permutation of the full f = {x,s) 
coordinates of any pair of particles - the \S^) spin state is required to be even with 
respect to the permutation of any pair of spin variables i and j. This condition limits the 
spin states to the manifold of maximal spin magnitude. Specifically, if we introduce the 
total spin operator, = Sj, then ■ [S*^) = y + l) |'S'^), corresponding 

to a total spin magnitude y, 
a linear combination of = y 



N ( N 
2 

i.e., all pseudo-spins aligned. Hence, the spin state is 
states of total spin projection Ms = — 

Y is an A^- 



Ms \S^ = Y,Ms). The M; 



s 



S 



N 



lt)l lt)2 - lt);v 



(76) 



\S^ = f,Ms) with S^\S^ = ^,Ms) 
spin stretched state. 

— Ms = — 
2 ' * 2 

In deriving the expression for the total energy E of the N interacting boson 
system, we use the Ising spin-spin form of the short-range particle-particle interaction. 
Integrating out the position variables explicitly, we encounter a volume v, the 'trap 
volume' that characterizes the 'tightness' of the confining Vp, 
1 



d X \xt [x) 



(77) 



For instance, the spin-independent interaction energy per particle, Cq, is inversely 
proportional to the trap volume v and proportional to the spin-independent interaction 
strength in the Ising form, Aq = {inh'^/m) 



4 



Ao 

Co = — • 

V 

The analogous spin-interaction energies per particle are given by the expressions 

AttH"^ / at — a\ \ 1 



(7J 



mo 

V 

2 (mil 



m 



2 



V m 
The total many-body energy, E, takes the form 



(a-^ + ai- 2au) 



(79) 



E = N 

-7 



-eo 



N 



N 



s 



N 



^^>+5^(5^|(-26^.,, + 6S,,)|^^> 



(80) 



so that the integration over the position variable with short-range interactions maps 
the N spin-1/2 boson problem into that of 1/2 spins coupled via an infinite range 
spin-spin interaction. By adding and subtracting 

^ AT 

^Z^^j,^^j> = (81) 
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we cast the Hamiltonian in terms of the total spin-operator = Ylj^=i^j- 
define effective magnetic fields that are c numbers, 

H, = (iV - 1) zeo 

HR = -2eHX + ez (82) 

The energy E, up to an unimportant shift, E' = E — N [ct + {N — 1) — then 
takes the form 



E'=(S 



■ (H, + H^) + 5f 5-;, (83) 

reminiscent of the Hamiltonian of magnetic single domain grains with anisotropic spin- 
spin interactions [38]. Chudnovsky and Gunther had pointed out that the anisotropy 
can set conditions under which we expect macroscopic quantum tunneling: sufficiently 
strong exchange interactions force the individual spins to align into a macroscopic spin 
vector, the anisotropy can give local energy minima corresponding to two distinct 
directions of the macroscopic spin and quantum mechanically, the macroscopic spin 
can travel through a classically forbidden region giving tunneling although the expected 
rate for such processes are exponentially suppressed with the number of spins. In the 
A^-boson quantum magnet, the alignment is enforced by permutation symmetry, the 
anisotropy caused by the Ising (or XY) nature of the effective inter-particle interactions 
and the number of bosons can, in principle, be controlled experimentally. 

We recognize the operator in the spin bracket of Eq. fl83l) as the spin Hamiltonian 
H. The Heisenberg equation of motion for the total spin operator, 

^-^^ ^S^ X [H„ + H^] + 26,Sf (S^ ■ z) , (84) 



ih- 



dt 

where Sj|^ denotes the part of the total spin vector that points in the direction of the 
magnetic field, S|^ = z (S^ ■ z), yields an Ehrenfest type of equation for the total spin 
expectation value if we take its expectation value. We can write the resulting equation 
as a Landau-Lifshitz equation (without damping term) [39], 

|(SN) = l(S^xHi„,,,) , (85) 

where the total magnetic field Utotai includes the Raman coupling effective field, the 
the effective, spin-independent short-range field carried by the other particles and the 
contribution caused by the Ising spin-spin interactions, 

Ht^tai = H„ + + e,Sf (S^ ■ z) . (86) 

In the absence of Raman coupling H/j = the total magnetic field Utotai points in 
the z-direction and all the spin expectation vector can do is precess around the z- 
direction. In fact, the derivative of the expectation value of any power of vanishes 
so that the conservation of up and down particles ensures that the distribution of the 
spin up and spin down particles will remain constant in time precluding collective 
tunneling of the spin in the absence of Raman coupling. A Raman pulse can then 
precisely control and initiate the macroscopic quantum tunneling while leaving the 
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other assumptions and parameters of the system untouched. In addition to tunnehng, 
the ground state of the system can be a superposition of two distinct states in each of 
which the total spin points in different directions. That Raman-coupled (or Josephson- 
coupled) two-component BEC systems can take on macroscopic Schrodinger cat states 
was pointed out in |10] and worked out in [H] - we simply determine the interaction 
parameters and indicate how these can be controlled. The Raman control provides an 
important advantage to the N-boson quantum magnet over the double-well proposals 
for realizing macroscopic quantum tunneling and creating macrocopic Schrodinger cat 
states. Varying the potential barrier in a double well system to control the tunneling 
can also render the two-state approximation invalid and lead to unwanted excitations. 

Observing the coherent oscillations of the total spin that is quantum tunneling can 
also test fundamental aspects of quantum mechanics (against macroscopic realism) by 
verifying Legget-Garg inequalities [12]. In addition, the A^-boson quantum magnet spin 
dynamics can also explore spin squeezing, non-classical quantum evolution near unstable 
trajectories [13] and, when the Ising interaction is eliminated by a Feshbach resonance, 
realize the Burnett-Holand proposal for Heisenberg limited interferometry [H] by using 
the Raman pulse as a beam-splitter. These connections become obvious using the spin- 
spin form of the inter-particle interactions which also reveal the control that routine 
cold atom knobs such as the intensity of the confining potential, the detuning of the 
Raman coupling pulse and the magnetic field of a Feshbach resonance can exercise. 

5. Conclusions 

In conclusion, we described the effective spin-dependent interactions of ultra-cold alkali 
atoms occupying two distinct hyperfine states in an external magnetic field. The 
magnetic field lifts the degeneracy of the atomic Zeeman levels and permits the selection 
of two hyperfine states to act as the effective 'spin-up' and 'spin-down' states of 
the particles so the atoms can mimic the behavior of magnetic spin-| particles. We 
described the spin-dependent effective interaction as a spin-spin interaction. The form 
of the effective spin-spin interaction depends explicitly on the quantum statistics of 
the interacting particles. As a consequence of the zero-range nature of the interaction, 
the interaction of spin-i bosons can be described as an Ising or, alternatively, as an 
XF-coupling. The parameters of the spin-spin interaction depend on the scattering 
lengths of the relevant binary alkali atom collision channels in the external magnetic 
field. For relatively low values of the magnetic field (sufficiently large to cause a 
Zeeman level splitting that permits the selection of two hyperfine levels) we calculated 
the parameters as a function of the external magnetic field in the Degenerate Internal 
State (DIS) approximation. We illustrated the advantage of the spin-spin interaction 
form by mapping the system of N spin-i bosons in a tight trapping potential on that 
of N spin-^ spins coupled via an infinite range interaction. The explicit expressions 
reveal which parameters of the spin Hamiltonian can be controlled and how. The spin 
Hamiltonian also suggests that the N-boson quantum magnet provides an intriguing 
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laboratory for the exploration of fundamental quantum studies. The list of promising 
uses include the study of collective quantum spin tunneling (which can be used for testing 
fundamental aspects of quantum mechanics), the controlled observation and utilization 
of spin squeezing and the creation and study of highly non-classical states. 
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